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Abstract. The iron-based superconductors allow for a zoo of possible order 
parameters due to their orbital degrees of freedom. These order parameters are often 
written in an orbital basis for a microscopic analysis and a comprehensive symmetry 
classification. Unlike in standard single-band superconductors, where electrons of the 
same band - and hence, energy - are paired, a general order parameter in such a 
multi-orbital system also contains pairing of electrons belonging to different bands. As 
this corresponds to pairing of electrons of different energy, such order parameters are 
energetically less stable. Here, we present a simple criterion for a stability analysis 
of the orbital part of the gap function based on the basic principle that electrons of 
equal energy are paired in the superconducting state. This not only allows to find 
the most stable states, but also to identify the terms in a (tight-binding) Hamiltonian 
suppressing any given superconducting state. Our approach thus allows to identify the 
minimal Hamiltonian necessary to compare competing instabilities. 



PACS numbers: 74.20.Rp, 74.70.Xa, , 74.62.Bf 



Gap Symmetry and Stability Analysis in the Multi-Orbital Fe-Based Superconductors 2 
1. Introduction 

The symmetry classification of superconducting order parameters lias proven to be 
a powerful tool in the context of unconventional superconductivity [1]. In single- 
band superconductors one distinguishes various gap functions A(k) according to their 
behavior under all the symmetry transformations of the generating point group of the 
crystal, i.e., the irreducible representation they belong to. While the interactions are 
the dominant factor for determining the pairing channel, many interactions allow for 
instabilities of different symmetry. With only gap functions belonging to the same 
irreducible representation mixing, this allows to find each class' critical temperature 
and thus the dominant pairing channel. An interaction with dominant nearest-neighbor 
character on a square lattice, such as antiferromagnetic spin fluctuations, for example, 
allows for extended-s-, or d-wave (singlet) order parameters, which can be analyzed 
independently for the leading inst ability |2]. 

In the Fe-based superconductors, several bands are involved in the low-energy 
physics. According to bandstructure calculations, the main contribution to the density 
of states at the Fermi level stems from Fe d orbitals hybridizing with As p orbitals[3],|4]. 
In most compounds, the bands crossing the Fermi energy have mainly dxz and dy^ 
orbital character and thus, minimal tight-binding models consider only these two 
orbitals[5],[6],[7]. After back-folding of the Brillouin zone due to the out-of-plane positions 
of the As ions leading to a two-Fe unit cell, these models indeed yield an overall 
agreement with the experimentally found Fermi surfaces. However, at least the three 
t2g orbitals have to be considered in order to find an agreement in the Fermi velocities 
as well [El [9]. 

This multi-orbital nature of the Fe-pnictides complicates any order parameter 
analysis: The additional orbital degrees of freedom lead to an abundance of gap 
functions, which are now matrices in orbital spac^. A symmetry classification of all 
possible gap functions with respect to the generating point group D^^ of the Fe-As layer 
has been done for two-, as well as three-orbital models for example in Refs. [TTl [T2] 
and Ref. [13], respectively. One difficulty arising with such an analysis is that the order 
parameter has to be written in orbital space, whereas the pairing of electrons in the BCS 
picture is best described in band space. Indeed, only pairing of electrons with the same 
energy, hence belonging to the same band, leads to a superconducting instability. While 
the analytical transformation into the band basis to analyze the stability of an order 
parameter is straight forward in a two-orbital description [TTl [T2| [T^. the generalization 
to three or even more orbitals is not possible. Alternatively, the loss of condensation 
energy due to nodes can be used as a criterion for the stability of the various gap 
functions. This requires knowledge of all the nodes of a given gap function and also 
only allows for an even cruder first estimate. 

In this work, we build on the aforementioned pairing of equal-energy electrons and 
present a simple criterion to analyze the stability of a gap with a given orbital structure. 

I This is similar to the spin-space matrices used when describing spin-triplet pairing [TO] 
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Simply speaking, we analyze whether a gap structure corresponds to a pure intra-band 
pairing, or whether electrons belonging to different bands are paired. For a systematic 
analysis, a comprehensive symmetry analysis of the free Hamiltonian and the order 
parameter, both matrices in orbital space, is necessary first. Our approach then also 
allows to identify the terms in the Hamiltonian that (partially) suppress the other gap 
structures. This has an important application in determining the minimal Hamiltonian 
needed to (numerically) compare different superconducting instabilities on equal footing. 
Note that such an analysis is unlike comparing extended-s- and (i-wave paring in the 
single-band case, where always electrons belonging to the same band are paired. Here, 
the suppression is due to inter-band pairing, more comparable to the suppression of 
spin-singlet superconductivity due to a magnetic field. 

After introducing the general scheme in the next section, we apply it to the case 
of the iron-based superconductors. In order to construct the general tight-binding 
models, we start in Section [3] with a few comments on the symmetry of the Fe-As 
layer and its generating point group D^h. We then focus our discusion in Section [Hon 
a two-orbital description of the Fe-As layer and, with experiments strongly supporting 
spin-singlet pairing, we restrict our analysis on this pairing channel p!5l [T6] . The two- 
orbital model carries the advantage that we can gain additional analytical insight from 
the linearized gap equation. Finally, we also comment on the three-orbital model in 
Section [5l As a starting point for both models, we present a comprehensive symmetry 
analysis of the Hamiltonian and the gap structures including the two-Fe unit cell. As 
an interesting result, we find that the extended-s and (ia.2_y2-wave pairing in the (two-Fe 
unit cell) three-orbital model is suppressed by the hybridization of the d^z/dy^ with the 
dxy orbitals, whereas this suppression is missing in general two-orbital models or when 
working with a one-Fe unit cell. 



2. Symmetry and Stability in Multi-Orbital Systems 

2.1. Single-Band Systems 

Before discussing the multi-orbital case, we give a brief overview over some symmetry 
aspects and notation for a superconductor in a single-band system. Only considering 
spin-singlet superconductivity, the mean-field Hamiltonian is a 2 x 2 matrix in Nambu 
space, = (Ckt^c^k;)' given by 

/MFn,^ _ ( ^S(k) A(k) 

A*(k) -Hi(-k) 

'Hj(k) = 'Ho(k) = .^k(= C-k) is the free hopping Hamiltonian|§| and A(k) = A?/'(k) is the 
superconducting order parameter or gap function. Since describes the non-interacting 
system, it has to transform trivially under time reversal and all the symmetry operations 
of the generating point group Q of the crystal. From the transformation behavior of 



§ For simplicity, we only consider time-reversal and inversion symmetric systems. 
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the single-particle operator under g & Q, gcy_ = c/j^k, and Rg the corresponding rotation 
matrix in momentum space, we find C,Rg]i = for all g E Q. In a similar way, the pair- 
wave functions 'i^(k) are symmetry classified in terms of the irreducible representations of 
Q, i.e., their behavior under ipCk) H- ip{Rg]<.). Note that the Pauli principle requires the 
full gap function to change sign under the exchange of the two electrons, leading to the 
well-known requirement that an even wave function ipCk) is combined with spin-singlet 
pairing. Within BCS theory, the mechanism of superconductivity requires the pairing of 
electrons of equal energy. In the above Hamiltonian, the gap function always connects 
two electrons with the same energy, thus allowing for a superconducting instability for 
arbitrary small (attractive) interactions. 



2.2. Multi-Orbital Systems 

In a system with n orbitals, the Hamiltonian and the (spin-singlet) gap function are x n 
matrices in orbital space, leading to a mean-field Hamiltonian analogous to Eq. ([T]) that 
is a 2n X 2n matrix, 

^MF.._/^0(k) A(k) \ 

^ ^"^^-U^k) -Ho(k)j' 

where we have again assumed time-reversal and inversion symmetry and for simplicity 
omitted any spin-orbit coupling. The Hamiltonian in orbital space is of the general form 

?i^o(k) = EA"/ia(k), (3) 

a 

where {A"} is a basis for the n x n matrices. The Hamiltonian Eq. ([2]) still has to 
transform trivially under all the symmetry operations of the generating point group Q. 
Since the symmetry transformations now act on the basis functions as gc]^ = UgCRg]f_, 
the individual elements of the Hamiltonian ([3]) transform as 

g[A^K{k)] = {U^A^Ul)K{Rgk). (4) 

For matrices and functions that transform as one- dimensional representations, simply 
the product of them has to transform trivially. For higher-dimensional representations, 
multiple matrices and momentum functions can be combined (Such an example is shown 
in Section [5]). It is thus helpful to write the Hamiltonian as 

?io(k)=<(k) + ^H^(k), (5) 

a 

where 7/Q°'(k) is proportional to the n x n identity matrix and "Hq (k) are all the other 
smallest parts transforming trivially. 

The (spin-singlet) superconducting order parameter A(k) can similarly be written 

as 

A(k)=5:A>„(k). (6) 
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This order parameter can now describe the pairing of electrons of all possible orbitals. 
As for the free Hamiltonian, the full symmetry of the order parameter can be found 
from 

^?[A(k)] = Y.{u;^^ul)URa^)■ (7) 

a 

The multi-orbital structure leads to a plethora of different gaps due to the additional 
degree of freedom. Most of these gaps will, however, lead to inter-band pairing and are 
thus less (or not at all) stable. After a brief discussion of the stability of general gap 
structures, we will see in Section H] for the two-orbital model and in Section [5] for the 
three-band model an explicit symmetry classification and stability analysis of such order 
parameters. 

2. 3. Stability of the Gap Structures 

In this section, we discuss a simple qualitative argument for the stability of gap 
structures and show that only the ones with ['Ho(k), A(k)] = are fully supported by 
the Hamiltonian. For a quantitative stability analysis of the terms in the Hamiltonian 
harming the superconducting state, one can for example resort to the linearized gap 
equation as done analytically for a two-band model in Section The superconducting 
instability requires to pair electrons of equal and opposite momentum, and most 
importantly, equal energy. For a gap structure that commutes with the orbital structure 
of the Hamiltonian, i.e. [Ho(k), A(k)] = 0, "Hoik) and A(k) can both be simultaneously 
diagonalized by the same n x n matrix U^. The block-diagonal matrix with the two 
blocks Uk thus leads to a diagonal form of the mean-field equation ([2]) with the gap only 
connecting states with equal energy, i.e., the gap function describes intra-band pairing 
only. If [?io(k), A(k)] ^ 0, we define the set C = {a| A(k)] = 0} and write the 
Hamiltonian as 

?io(k) = nifik) + y: K(k) + y: ^^(k) (8) 

with the first sum commuting with A(k), while the second sum does not. If the 
Hamiltonian consisted of only the first two terms, the gap structure would describe 
intra-band pairing only and would be fully supported. The last sum, however, leads to 
inter-band contributions. The fact that there are terms in this last sum does, however, 
not necessarily lead to a full suppression of the instability. As we will see explicitly for 
the two-band model, the stability instead depends on the respective energy scales of 
the two groups of terms. This scheme thus allows for a first estimate of the stability 
of all possible gap functions and also provides an insight into what terms in the non- 
interacting Hamiltonian harm a specific superconducting gap function. 

3. Symmetry of the Fe-As Layer 

In order to apply the above scheme, we first recapitulate the symmetry properties of a 
single Fe-As plane [HI [I3]. Figure [I](a) shows the structure of the Fe-As layer with its two- 
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I 



Figure 1. (a) Cristal structure of the Fe-As layer with the two Fe sublattices A (dark 
circles) and B (white circles) . The small crosses (dots) in the middle of the Fe squares 
are the As sites above (below) the plane. Also, the symmetry operations belonging 
to D2d are depicted with C2 and ^4 being around the square on one of the Fe sites. 
Finally, the star between two iron sites denotes a center of inversion, (b) Orbital basis 
for the two-Fe unit cell. 



sublattice structure due to the out-of-plane positions of the As ions. The individual (two- 
Fe) unit cell is invariant under the following symmetry operations: The identity E, a 
two-fold rotation C2 and two improper rotation 5*4, both around the z axis, two two-fold 
rotations around the x (y) axis, as well as two mirror planes ad along the diagonals 
of the xy plane, together building the point group = {£", C2, 2C2, 25*4, 2(7^}. Note 
that the generating point group of the layer is isomorphic to -D4/1, where all the 
elements 1^4/1 \ D2d have to be combined with a lattice translation interchanging the 
two sublatticedjll- Taking this into account ensures having the inversion as a symmetry 
element of the layer, though not the unit cell. This is crucial for both, the construction 
of the Hamiltonian as well as the classification and analysis of the superconducting order 
parameter. 



4. Two-Orbital Model 

Hopping Hamiltonian 

In this section, we discuss the symmetry allowed terms of the multi-orbital hopping 
Hamiltonian based on the d^z and dyz orbitals with real-space operators d^^/^ and d^^f\ 
respectively. We work with the two-Fe unit cell, and hence, additionally introduce the 
states d!fg^^ and d^is^\ where the orbitals on the B sublattice are multiplied by —1 (see 
Fig. [Hb)). In momentum space, this translates to having two operators d^^^ / (ik+Q^ 
and d^g^ / c^k+Qs; respectively with Q = (7r,7r). It is thus convenient to work in a 
tensor space of the {d^z, dyz) orbital space and the (k, k -|- Q) momentum space. The 

II A crystal with this property is called non-symmorphic. 
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IR(A/.) 


^19 




Big 


B2g 


Am 


A2u 


Blu 


B2u 


Y^ab 



















Table 1. Classification of the orbital basis matrices S"'' ~ p"' ® t'^ for a two-orbital 
dxz/dyz basis. Note that the behavior under inversion (u/g) solely depends on the 
momentum part (r^), while the behavior under the symmetry transformations of the 
unit cell {D2d) is given by the orbital part (see [Appendix K\ . 



IR 




f3 




1, COS kx cos ky, cos 2kx + cos 2ky 


COS A;^: + cos ky 


Big 


cos 2kx — cos 2ky 


cos /c^ — cos ky 


B2g 


sin kx sin ky 










Eu 


{sin kx cos ky, sin ky cos /Ca;} 
{sin2/i;a;, sm2ky} 


{sin /cj., sin ky} 



Table 2. The basis functions for hopping terms /lab(k) up to third-nearest neighbors 
classified in terms of irreducible representations of D/^h- The functions are divided into 
the ones that change sign for k i— )■ k -I- Q, thus allowing for combinations with f'^ (f^) 
in the Hamiltonian, and those that do not, thus allowing for combinations with f° 
(fi). 



Hamiltonian can then be factorized in analogy to Eq. ([3]) as 

^o = Ecl[ES'^''^«''(k)]c;, (9) 

k,s a,b 

with Cks = (c^kT^ d^Vs\ '^k+Qs5 ^k+Qs)^ ^^iid S"'' = p"-®f^ is a tensor product of and f'', 
the identity (a, fe = 0) and Pauli matrices {a,b = 1,2, 3) for the {d^z, dyz) and (k, k -|- Q) 
spaces, respectively. In the following, we separately analyze the behavior of /iab(k) and 
S"'' under all the symmetry transformations of the Fe-As layer. 

Table [1] shows the classification of all the matrices S"'' with respect to D^h (see 



Appendix A for details). As can be seen from this table, the 'momentum' part f 



b 



determines the behavior under inversion (u/g). This has important consequences for 
the possible terms that can appear in the Hamiltonian: Momentum functions that are 
even (odd) under inversion can only be combined with f'' and (f^ and f^). Table 
|2] summarizes all the hopping factors /ia6(k) from up to third- nearest-neighbor hopping 
together with the f matrices in momentum space they have to be combined with. By 
construction, also the transformation k i— )■ k -|- Q has to be considered, such that for 
example the basis functions (cosfc^: ± cos ky) = —[cos{kx + ir) i: cos{ky + vr)] have to be 
combined with f^. 

We are now in the position to construct the Hamiltonian by requiring it to transform 
as the trivial representation Aig. Since only one-dimensional representations appear in 
Table [1], two elements of the same irreducible representation have to be combined, e.g. 
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s-wave 



(b) extended- 



(d) 



(e) 
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E*^" E'''^(cos fc^ + cos E"" cos fc^^ cos fcy E^^(cosA;j: — cosfcj^) E^° sinfcj; sinfej, 



Figure 2. Sketch of the possible spin-singlet gap structures in momentum space with 
the full Aig symmetry, (a) shows the standard s-, (b) extended-s- and (c) s±-wave 
gap functions, (d) and (e) are possible fully symmetric gap functions with band- and 
momentum structure belonging to Big and i?2g, respectively. The two colors depict 
the sign change of the gap and the dashed (dotted) lines the nodes of the momentum 
(orbital) part of the gap. Note that for the gap functions belonging to Aig these two 
lines have to fall together. 



COS COS ky with The Hamiltonian reads 

no = Ecl[S°°/^oo(k)+S°^/io3(k)+S^°/i3o(k)+S^^/i33(k)+S^°/iio(k)]Ck,, (10) 

k,s 

with /loo(k) = 4:t2s cos kx cos ky + 2t3s(cos 2/^^. -|- cos2ky) — /i, ho^ik) = 2tis(cos/c^ + 
cos ky), /i3o(k) = 2t3ii{cos2kx — cos2ky), /i33(k) = 2tid{coskx — cos ky), and /iio(k) = 
4:t2d sin kx sin ky . Note that these hopping terms correspond to isotropic hopping along 
the X and y axis (/loo(k) and /io3(k)), anisotropic hopping with respect to the axes 
(/i3o(k) and /i33(k)), as well as inter-orbital hopping (/iio(k)). The Hamiltonian (fTOl) has 
4 bands with dispersions 

d"^(k) = ei\k) ± v/[eL"HkP + [/^io(k)]2, (11) 

with e^+^^\k) = /loo(k) ± /io3(k), e'L^^\k) = h^oik) ± h-ssik), and a = (1,2) for the 
subspace of (k, k -|- Q). Note that no terms transforming trivially can be constructed 
containing the matrices or and thus, there is no allowed coupling terms between 
k and k + Q. The first Brillouin zone is thus only artificially folded down to the two-Fe 
zone. Terms that explicitly fold the Brillouin zone in such a two-orbital description 
necessarily involve spin-orbit coupling. When sketching the gap functions below, we 
thus consider only one of the two subspaces for simplicity. 



4-2. Gap Glassification and Stability 

Only considering spin-singlet gap functions already allows for a large variety of different 
gap functions (see e.g. [HI [I2]). The combination of all the matrices from Table [H 
with momentum functions from Table [2] allows for 54 different order parameters and 
we thus discuss only some examples. Figure |2] sketches the Aig gap functions with 
momentum parts going up to next-nearest neighbors and only showing one of the 
momentum subspaces in the one-Fe Brillouin zone. In addition to the conventional 
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(c) 




E (cos kx — cos k, 



S"" sin k^ sin k,. 



(d) 
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I J 1 1, 

W ! ^ 
















COS kx COS fcj, 



Figure 3. Momentum space representations of the gap functions with Big [(a), (b), 
and (c)] and B^g [(d) and (e)] symmetry, (a) is the standard dx2-y2- and (d) the dxy 
wave gap. (b), (c), and (e) have an Aig momentum dependence combined with a Big 
and B2g band dependence, respectively. In addition to the nodes dictated by their 
momentum part (dashed hues), they thus also have nodes coming from their orbital 
structure (dotted lines). 



s-wave, extended-s and s± gap functions, given by 

A(k) = (12) 



A(k) = S°=^(cos K + cos ky), (13) 

A(k) = ±^^ cosh cos ky, (14) 
also gap functions with gap nodes on the Fermi surface around the F point, namely 

A{k) = t^^ {cos k:,- cos ky), (15) 

A{k) = t^^ sin k^ sin ky (16) 

are shown. 

Figure [3] shows gap functions that transform as and i?2c,- These are, first, the 
so-called dx^-y^-wave gap (Fig. Et^a)) 

A(k) = S°2(cos kx - cos ky) (17) 

and dxy-wave gap (Fig. EJj^d)) 

A{k) = t^^ sin kx sin ky. (18) 

In addition, three gap structures are shown with a 'trivial' momentum part in 
combination with a B^ (Fig. |3]^b), (c)) and i?2g (Fig. EJ^e)) orbital part, respectively, 

A(k) = S^2(cos kx + cos ky), (19) 
A(k) = cos kx cos ky, (20) 
A(k) = S^°cosA;^.cosA;y. (21) 

The gap structures with non-trivial orbital part in Eqs. ( |T5l) and ( fT6ll . as well as (fT9ll 
- fl2T]) have not only gaps due to their momentum dependence (dashed lines), but also 
their orbital dependence (dotted lines). Note, however, that here the former ones are 
not protected by symmetry and thus any admixture of other gap structures belonging 
to the same irreducible representation can shift or open these nodes. Nevertheless, these 
gaps by themselves are not fully gapped and hence, in general less stable. 
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The fact that the gaps with non-trivial orbital structure have (additional) nodes 
and thus a reduced superconducting condensation energy already points toward less 
stable superconductivity. Even without knowledge of these gap nodes, we can analyze 
the stability of the above gap structures using the criterion introduced in Section [2j 
All the gap structures that transform trivially in the two subspaces (k, k + Q), namely 
the ones defined in Eqs. f|T2|) - f|T4|) . f|T7|) . and ( !T8|) . commute with the free Hamiltonian 
ffTOj) . These gap structures are thus intra-band only and fully supported. However, 
the gap structures with non-trivial orbital parts do not commute with the Hamiltonian 
(fTOj) . The gap structures defined in Eqs. ( !T5|) . ( IT9|) . and (!20|) do not commute with the 
term S^'^/iio(k) in the Hamiltonian (11 01) . while they commute with all the other terms. 
This means that these gap functions are suppressed by inter-orbital hopping. The gap 
structures in Eqs. ( fT6l) and ( l2Tl) . on the other hand, do not commute with the terms 
S^'^/i3o(k) and E^'^/i33(k), i.e., they are suppressed by the 'one-dimensional nature' of 
the I dyz bands. 

Finally, we comment on the peculiar gap structures that combine an odd momentum 
function with spin-singlet pairing. In order to obey the Pauli principle, either or 
has to appear in the decomposition of S"^. Examples are 

A(k) = S^^sinfc^, (22) 

belonging to the two-dimensional representation Eg, or 

A(k) = S^^sinfc^ (23) 

belonging to Eu- The former gap structure corresponds to finite- momentum pairing in 
the one-Fe zone and it is thus suppressed. The latter gap structure only commutes with 
the terms proportional to S*^'' and in the Hamiltonian ( ITOl) . From our discussion so 
far, we hence also expect this gap structure to be suppressed. Calculating the linearized 
self-consistency equation in the next section, we can explicitly analyze how the various 
terms influence the pairing instability for the above gap structures. 



4-3. Linearized Gap Equation 

In this section, we evaluate the linearized self-consistency equation for the gap structures 
to see explicitly how the criterion introduced in Section |2] works out. For this purpose, 
we introduce a spin-singlet pairing interaction of the general form 

= T7 X! ^a/3,Aii^(k, k )C^j^^C^_i^|C^_j^/^C^l^/^ (24) 

k,k' 

(We use the convention of summing over repeated indices). In the framework of the 
Gor'kov equations, we find the linearized gap equation 

A^pik) = -T5^5^V,/3,M.(k,k')x[Go(k',u;„)A(k')G;^(-k',-a;„)] ^,(25) 

where (jo(k, cj„) = [iw^S^^ — ?^o(k)]~^ is the non-interacting Green's function with 
Un = (2?7, + l)7rT the fermionic Matsubara frequencies. This equation determines the 



Gap Symmetry and Stability Analysis in the Multi- Orbital Fe-Based Superconductorsll 



critical temperature T^, for a given gap structure, where the largest belongs to the 
leading instability. Parametrizing the interaction similar to Eq. ([6]), 



a,b I 



\4^,,.(k,k') = EE^'i'FVi;^(k)J^^ X [S'^Vll^(k) ,„ , (26) 



with ip'^l (k) the individual basis functions from Table [21 we find the self-consistency 
equation for each individual gap structure. 

We neglect here all possible gap mixing and only focus on the diagonal elements. 
We first consider gap structures that commute with the Hamiltonian, namely the ones 
with a 'trivial' structure S"*^ and S"^, respectively, given in Eqs. flT2|) - f|T^ . ffTTj) . and 
( |T8|) . After summation over the Matsubara frequencies, the linearized gap equation for 
each gap reads 



where ah = 00,03 and ip'^l{\<.) is the momentum part of the respective gap. In this 
equation, every state with energy close to the Fermi energy adds to the instability due 
to the intra-band (equal-energy) character of the pairing with an instability even for an 
infinitesimal attractive interaction due to the divergence for T — )■ 0. As expected, these 
gap structures are thus fully supported by the hopping Hamiltonian. 

The non-trivial gap structures from Eqs. (fT5|) and (fT9|) (and analogously for 
Eq. (l20l) ) yield a self-consistency equation of the form 



'\hJ[hioik'W + [e^-W 2e«(k') '""H 2T 



[hoW [V'g(k')]\„„wC(k' 



tanh(S^)}. (2, 



Note that the first term on the right-hand side again corresponds to intra-band pairing 
with the T — )■ divergence, whereas the second term comes from inter-band pairing 
and has no such divergence. Hence, this self-consistency equation still allows for a 
solution for an infinitesimal interaction, however, Tc is suppressed. As expected from 
the discussion in the previous section, the ratio hio(k) / e"L(k) determines the strength 
of the instability. Analogously, we find 



[/.lo(k')]^ [<^(k')]\„„wC(k' 



— V 



(0 

'\f:^^^[ho{k'W + [eU^')? 2e(k') ^^""V 2T 



^ I \h.JV'Wi ^ rkM12 2P^{V'\ 



tanh(i^)} (29) 



, [e4kO]^ [^S(kO]^ , , /Ok' 

[/iio(k')]2 + [£^(kO]2 24(kO V 2T 

for the gap functions from Eqs. ( fT6l) and ( 12T]) with the same structure as Eq. ([28 
Finally, for the gap structure given in Eq. (!23l) we find 



23 



This gap function is only describing inter-band pairing and thus has no T — t- divergence 
and associated instability for an arbitrarily small interaction any more. 
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AO, A8 


A2 


A^ 


Bi 


A=^ 


B2 


Ai 


E 


{A^A6}, {A^A^} 



Table 3. Classification of the Gcll-Mann matrices in terms of irreducible 
representations of D2d- The momentum part then determines how the matrices 
A"'' = Aq ® fb behave under inversion and thus what irreducible representation they 
belong to (see [Appendix A| for details). 



5. Three-Orbital Model 



We now turn to the three-orbital model, where in addition to the d^z and dyz orbitals, 
also the d^y orbital is introduced. This is necessary for an accurate description of the 
low-energy electronic structure and we hence briefly repeat the analysis of the symmetry 
and stability of the gap functions. The Hamiltonian can again be factorized as 

^o = ECl[EA"'^a.(k)]4, (31) 

k,s a,b 

where now, 

r _ (A^^^ ri^y^^ ri^^y^ ri^y^^ A^y^ r^o^ 

L-ks — l"ks ' "ks ; "ks ' "k+Qs; "k+Qs' "k+QsJ • l-J^J 

We have introduced the matrices A*^^ = A*^ ® f'', with as before the matrices in 
{k, k+ Q} space and the orbital basis is given in terms of the (3 x 3) Gell-Mann matrices 
A" (see Appendix B ). In analogy to Section UT] we analyze the behavior of these matrices 
under the symmetry operations of the generating point group. For simplicity. Table [3] 
only lists the classification of the Gell-Mann matrices. The full matrix A"'' = A'^®f'' then 
transforms like A'^ under and is even (odd) under inversion for 6 = 0,3 {b = 1,2), 
c.f.. Tabled] All the matrix elements only connecting the dxz and dyz orbitals are (with 
the appropriate matrices) the same as for the two-orbital case. For the diagonal part of 
dxy, we simply add the terms 

i(A«»-V3A«)2i.,(cos*. + cos*.) (33) 

+ ilA" - V3A*°)(4«;„cos«:,cos*, - (34) 

including nearest- and next-nearest-neighbor hopping and the chemical potential fi^y 
Due to the out-of-plane position of the As atoms, there is also an allowed hopping 
between the dxz / dyz and the dxy orbitals. This is a nearest-neighbor hopping of the 
form 

A^2 sin kx - A^2 sin ky (35) 
and a next-nearest-neighbor hopping 

A^"^ sin kx cos ky — A''"^ sin ky cos kx, (36) 
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both odd under inversion. Note that these terms now are made from matrices / 
momentum functions that transform hke the two-dimensional representation and 
can thus not be spht further. To choose which Eu matrices from Table [3] have 
to be used, time-reversal symmetry has to be considered, which requires HQ^Ck) = 
[iJo^(— k)]*. Further note that the Hamiltonian again separates into two blocks, with 

(4T\4T\4TqJ and vice versfl 

With this hopping Hamiltonian, we can again analyze the stability of possible gap 
structures. Compared with the two-orbital case, there is one key difference: Even 
without spin-orbit coupling the k and k -|- Q subspaces are connected through the 
hopping terms in Eqs. ( l35l) and ( 1361) . As a result, the only spin-singlet gap function that 
has no inter-band contributions and hence is fully supported by the Hamiltonian has 
the trivial orbital structure A°°. This has two important consequences: First, the most 
stable gap structure has the same momentum dependence on all Fermi surfaces, i.e., 
there is a single gap function determining the gap everywhere in the Brillouin zone [13] • 
Second, the orbital structure only allows for functions '0(k) that transform trivially 
under k i— t- k + Q. This means that the gap function with extended-s and dx'^-y'^ 
symmetry with = coskx ± cos ky are suppressed by the nearest-neighbor hopping 
between the dxz / dyz and the d^y orbit alsE] 

6. Discussion and Conclusion 

In this work, we have introduced a simple criterion to analyze the stability of the 
orbital structure of superconducting order parameters in the context of multi-orbital 
systems: A vanishing commutator of the gap function with the free Hamiltonian. This 
argument is based on intra-band pairing, i.e, the equal-energy pairing of electrons 
in a weak-coupling approach. Depending on the structure of the interaction, the 
instability can therefore still occur in an orbital channel that is not fully supported by 
the Hamiltonian, especially with almost degenerate bands and within strong coupling. 
The analysis is thus only a first step and resembles the argument regarding the nodes 
in the superconducting gap due to non-trivial orbital structures and the associated 
loss of condensation energy. However, the introduced scheme additionally allows for 
the identification of the superconductivity-suppressing hopping terms. In the case 
of competing gap functions, this can help finding a minimal model for an unbiased 
description. 

We have applied our scheme to the case of spin-singlet pairing in the Fe-based 
superconductors. Looking at all symmetry-allowed hopping terms in a three-orbital 
model for the two-Fe unit cell, we have found that only the trivial orbital structure 
proportional to the identity is fully supported by the system, i.e., only describes intra- 

% This separation into two blocks can be understood from the behavior under a glide-plane 
transformation as noted by Ref. [5] 

+ Note that this suppression is present even without hybridization of the electron pockets discussed in 
Ref. [n]. 
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band pairing. In particular, the gap functions with ip(k+ Q) = —ipik) are suppressed 
by the hopping terms connecting k and k + Q. Moreover, the gap functions belonging to 
Big in a two-dimensional model have additional gap nodes due to their orbital structure 
along the folded Brillouin zone, as this symmetry belongs to B2g in this reduced Brillouin 
zone[T8]. The gap function shown in Fig. [3](c) then corresponds to a 'd±' gap[T9]. to 
give an example. As the suppression and the additional nodes are missed considering 
two-orbital or general one-Fe-unit-cell models, our analysis emphasizes the importance 
of using at least a (two-Fe) three-orbital description when studying possible competing 
order parameters. Especially in the description of Fe-based superconductors with only 
electron pockets this could be relevant [201 |2T[ 122] . 

Finally, we comment on the generalization of the introduced scheme to spin- 
triplet superconductivity and the introduction of spin-orbit coupling. For simplicity, 
we have focused in this work on spin-singlet pairing and no spin-orbit coupling. For a 
generalization, the Nambu spinors have to be written in a way as to preserve the form 
of the mean-field Hamiltonian ([T]). The application of the commutation criterion is then 
again straight forward. In the case of a single-band system with time-reversal symmetry, 
for example, the appropriate basis is ^(k) = (ckt, Ck|, c^kj,' The lifting of the 

degeneracy of the spin-triplet gap functions in non-centrosymmetric superconductors [23] 
is for example easily obtained in this way. 
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Appendix A. Transformation of the Basis States 

In this appendix, we explicitly analyze how the symmetry transformations of the 
generating point group of the Fe-As layer, D^h, act on the basis states 

of the two-orbital model introduced in Section |H For this purpose, we first analyze 
the behavior under all the symmetry operations of the unit cell D2d and than add 
inversion, since D^^ = B>2d ® I- It is convenient to look at the respective states in 
real space, as depicted in Figure Ht^b). The two-fold rotation around the z axis, C2, 
amounts to multiplication of all the states by —1, i.e., the transformation matrix is 
given by = — S^". Therefore, all the basis matrices transform trivially under C2, 
i.e. JJc2^°'^Uc2 = This already means that all possible orbital matrices only belong 
to one-dimensional representations of D2d, as can be seen from the character table lAll 

(xz) 

We continue with the two-fold rotation around the y axis. Under such a rotation, dig 
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IR 
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2S4 


C2 


2C'2 


2ad 
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2C4 




2(T^ 


2C'^ 




1 
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1 
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1 


1 


1 


1 
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^9 


1 


1 


1 


-1 


-1 


1 


1 


1 


-1 
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-1 


1 
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-1 


1 


-1 


1 


1 


-1 
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-1 


1 


-1 


1 


1 


-1 


1 


-1 


1 


E9 


2 





-1 








2 





-1 
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-1 

-1 
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-1 
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-1 

-1 


-1 

-1 


-1 

1 


-1 

-1 


-1 

-1 


-1 

1 
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-1 


1 


-1 


1 


-1 


1 


-1 


1 


-1 


Blu 


1 


1 


1 


1 


1 


-1 


-1 


-1 


-1 


-1 


B2U 


1 


1 


1 


-1 


-1 


-1 


-1 


-1 


1 


1 


Eu 


2 





-1 








-2 





1 









Table Al. Character table for 1)4/1 = ® I- Note that the symmetry operatfons 
are ordered such that the first five columns belong to £>2d. 

{d\g^^) stays the same, while d[^^^ {d[l^^) is multiphed by —1. The improper rotation ^4 
is a combination of a C4 rotation around z, followed by a mirror operation cr^ at the xy 
plane. This will take df^^^ to d'l^^\ d^^^"* to —dlg^\ as well as df^^^ to —di^^^ and d^l^^ to 
-djf \ Finally, the mirror operation can be obtained by combining C'2 and S4. This 
then yields the following eight matrix representations. 

By calculating how the matrices tl""^ = Pa® T'b transform under these transformations, 
we can thus analyze what irreducible representation of D2d they belong to. 

Finally, the Fe-As lattice possesses centers of inversion lying between neighboring 
Fe sites, which transform the basis as f// = TP^ . The transformation properties of the 
individual basis matrices under D2d and inversion leads to the classification shown in 
Table [D Note that S'^^ with b = 1,2 are odd under inversion. As can be seen from 
Table IXTl this means that for example which belongs to Ai in D2di belongs to Biu 
in -D4/1. This can be understood from the transformation behavior of the f matrices |23]: 
While f° and transform trivially (^ig), the matrices and transform as Biu in 
D4/1. As shown in Table \K\\ this is the irreducible representation that transforms in the 
same way as Aig in D2d- The total matrix S"'' then transforms as the product of the 
irreducible representation of the orbital part with Bi^- 

Appendix B. Gell-Mann Matrices 

For the three-orbital model used in Section |5l we introduce the Gell-Mann matrices: 
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.(B.l) 



Note that the matrices (2Ao + VSXs)/'^, Ai, A2, and A3 correspond to the (2 x 2) identity 
and Pauh matrices of the two-band model and 

/ \ 



-(Ao - a/SAs 




\ 1 



(B.2) 
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